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COMPUTING  RADIATIVE  HEAT  TRANSFER 
OCCURRING  IN  A ZONE  FIRE  MODEL 


Glenn  P.  Forney* 


Abstract 

Radiation,  convection  and  conduction  are  the  three  mechanisms  which  a zone  fire 
model  must  consider  when  calculating  the  heat  transfer  between  fires,  wall  surfaces  and 
room  gases.  Radiation  dominates  the  other  two  modes  of  heat  transfer  in  rooms  where 
there  are  fires  or  hot  smoke  layers.  The  computational  requirements  of  a radiation 
model  can  also  easily  dominate  the  work  required  to  Ccdculate  other  physical  sub- 
models in  a zone  fire  model. 

This  report  presents  algorithms  for  efficiently  computing  the  radiative  heat  ex- 
change between  four-wall  surfaces,  several  fires  and  two  interior  gases.  A two-wall  and 
a ten-wall  radiation  model  are  also  discussed.  The  structure  of  this  radiation  model  is 
exploited  to  show  that  only  a few  configuration  factors  need  to  be  calculated  directly 
(two  rather  than  16  for  the  four-wall  model  and  eight  rather  than  100  for  the  ten-wall 
model)  and  matrices  needed  to  solve  for  the  net  radiative  flux  striking  each  surface 
are  shown,  after  the  appropriate  transformation  is  taken,  to  be  diagonally  dominant. 
Iterative  methods  may  then  be  used  to  solve  the  linear  equations  more  efficiently  than 
direct  methods  such  as  Gaussian  elimination. 

The  radiation  exchange  algorithms  are  implemented  as  FORTRAN  subroutines 
named  RAD2,  RAD4  and  RADIO.  These  subroutines  along  with  a test  driver  are 
available  from  the  author. 


1 Introduction 

1.1  Background 

Radiation  is  an  important  heat  transfer  mode  to  represent  in  a zone  fire  model  due  to  the 
high  temperatures  attained  in  rooms  with  fires  or  hot  smoke  layers.  It  can  easily  dominate 
convective  and  conductive  heat  transfer.  A radiative  heat  transfer  calculation  can  also  easily 
dominate  the  computation  in  any  fire  model.  This  is  because  radiation  exchange  is  a global 
phenomena.  Each  portion  of  an  enclosure  interacts  radiatively  with  every  other  portion  that 
it  ‘sees’.  It  is  therefore  important  to  design  radiation  exchange  algorithms  that  are  efficient 
as  well  as  accurate. 

Most  zone  fire  models  use  two  wall  segments  to  model  radiation  exchange.  Harvard 
V [1],  FIRST  [2],  BRI  [3,  4]  and  FAST  [5]  are  some  examples.  FAST/FFM  [6]  on  the 
other  hand  uses  many  surface  segments  in  order  to  model  the  radiative  interaction  between 

‘Building  cind  Fire  Research  Laboratory,  Nationeil  Institute  of  Standards  eind  Technology,  Gaithersburg, 
MD  20899,  U.S. A.  (gpfrn8cfr6.cfr.nist.gov). 
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wall  surfaces  and  furniture  elements.  Harvard  V and  FIRST  model  the  two  wall  segments 
as  infinite  parallel  plates.  BRI  and  FAST  model  the  two  wall  segments  as  an  extended 
floor  and  ceiling.  The  extended  ceiling  consists  of  the  ceiling  plus  the  four  upper  walls^ 
Similarly,  the  extended  floor  consists  of  the  floor  plus  the  four  lower  walls.  The  purpose  of 
the  work  described  in  this  report  then  is  to  enhance  two  wall  radiation  exchange  algorithms 
by  considering  more  wall  segments.  In  particular  for  the  four-wall  case,  this  allows  the 
ceiling,  the  upper  wall  segment,  the  lower  wall  segment  and  the  floor  to  transfer  radiant 
heat  independently. 

This  report  describes  three  algorithms  for  computing  radiative  heat  transfer  between 
the  bounding  surfaces  of  a compartment  containing  upper  and  lower  layer  gases  and  point 
source  fires.  The  first  algorithm  uses  two  wall  segments,  the  second  uses  four  wall  segments 
and  the  third  uses  ten  wall  segments.  These  algorithms  each  use  the  net  radiation  equation 
as  described  in  Siegel  and  Howell  [7,  Chapter  17]  which  is  based  on  Hottel’s  work  in  [8]. 
An  enclosure  is  modeled  with  N wall  segments  (for  our  case  N will  be  2,  4 or  10)  and  an 
interior  gas.  A two  layer  zone  fire  model,  however,  requires  treatment  of  an  enclosure  with 
two  uniform  gases  (the  upper  and  lower  layer).  Hottel  and  Cohen  [9]  developed  a method 
to  handle  this  case  by  dividing  an  enclosure  into  a number  of  wall  and  gas  volume  elements. 
Treatment  of  the  fire  and  the  interaction  of  the  fire  and  gas  layers  with  the  walls  is  based 
upon  the  work  of  Yamada  and  Cooper  [10,  11]  on  N-wall  radiation  exchange  models.  They 
model  the  fire  as  a point  source  of  heat  radiating  uniformly  in  all  directions  and  use  the 
Lambert-Beer  law  to  model  the  interaction  between  heat  emitting  objects^  and  gas  layers. 

The  number  of  wall  segments  influences  the  accuracy  and  the  efficiency  of  the  radiation 
exchange  model.  The  suitability  of  a two,  four  or  ten  wall  radiation  model  then  depends 
on  the  particular  application.  More  wall  segments  are  needed  to  model  radiation  exchange 
accurately  in  rooms  with  large  temperature  variations. 

The  two,  four  and  ten  wall  algorithms  are  implemented  as  FORTRAN  77  subroutines 
named  RAD2,  RAD4  and  RADIO.  The  routines,  RAD2  and  RADIO,  take  advantage  of  the 
modular  structure  of  RAD4  by  using  a number  of  its  routines.  It  should  be  pointed  out 
that  the  computational  requirements  of  a general  N-wall  radiation  model  are  too  great  for 
now  to  justify  incorporating  it  into  a zone  fire  model.  By  implementing  the  net  radiation 
equation  for  particular  N (two,  four  or  ten  walls),  significant  algorithmic  speed  increases 
were  achieved  by  exploiting  the  structure  of  the  simpler  problems. 

1.2  Overview 

To  put  the  two,  four  and  ten  wall  radiation  models  into  perspective,  an  N-wall  model  is 
presented  using  the  notation  found  in  [7,  Chapter  17].  Modeling  assumptions  are  given  and 
the  net  radiation  equation  is  derived  for  the  general  N-wall  segment  case.  The  N-wall  model 
is  then  adapted  for  two,  four  and  ten  walls.  Two  key  results  allow  for  significant  algorithmic 
speed  ups.  First,  it  is  shown  that  only  eight  configuration  factors  need  to  be  calculated 
directly  rather  than  100  for  the  ten-wall  case  and  two  configuration  factors  rather  than  16 
for  the  four-wall  case.  This  is  important  because  the  configuration  factor  calculation  is  one 
of  the  major  bottlenecks  in  computing  radiation  heat  exchange.  Second,  it  is  shown  how  a 
matrix  involved  in  the  solution  of  a set  of  simultaneous  linear  equations  can  be  transformed 
into  a diagonally  dominant  matrix.  This  is  significant  because  an  o(A^)  rather  than  an 

^The  upper  wall  is  the  portion  of  a wall  above  the  layer  interface.  Likewise,  the  lower  wall  is  below  the 
interfeice. 

^fires,  walls  or  gas  layers  for  excimple 
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o{N^)  algorithm  can  then  be  used  to  solve  the  resulting  linear  system.  This  savings  is  more 
important  for  the  ten-wall  case  (or  larger)  than  the  four-wall  case.  Some  computational 
results  are  presented  and  the  subroutines  RAD2,  RAD4  and  RADIO  are  documented.  The 
structure  of  the  subroutines  used  with  RAD4  is  documented  and  a short  description  of  each 
subroutine  used  by  RAD2,  RAD4  and  RADIO  is  given. 

2 The  Problem 

2.1  N Wall  Segment  Radiation  Exchange 

The  N-wall  radiation  model  described  in  this  section  considers  radiative  heat  transfer  be- 
tween wall  segments,  point-source  fires  and  two  gas  layers.  An  enclosure  or  room  is  par- 
titioned into  N wall  segments  where  each  wall  segment  emits,  reflects  and  absorbs  radiant 
energy.  The  interior  of  the  enclosure  is  partitioned  into  two  volume  elements;  an  upper  and 
a lower  layer.  The  problem  then  is  to  determine  the  net  radiation  flux  emitted  by  each 
wall  segment  and  the  energy  absorbed  by  each  layer  given  the  temperature  and  emittance 
of  each  wall  segment  and  the  temperature  and  absorptance  of  the  two  gas  layers. 

These  calculations  can  be  performed  in  conjunction  with  a zone  fire  model  such  as 
CFAST.  Typically,  the  solution  (wall  temperatures,  gas  layer  temperatures  eic)  is  known  at 
a given  time  t.  The  solution  is  then  advanced  to  a new  time,  t + At.  The  calculated  radiation 
fluxes  along  with  convective  fluxes  are  used  as  a boundary  condition  for  an  associated  heat 
conduction  problem  in  order  to  calculate  wall  temperatures.  Gas  layer  energy  absorption 
due  to  radiation  contributes  to  the  energy  source  terms  of  the  associated  zone  fire  modeling 
differential  equations.  The  time  step.  At,  must  be  chosen  sufficiently  small  so  that  changes 
in  wall  temperatures  are  small  over  the  duration  of  the  time  step. 

2.1.1  Modeling  Assumptions 

The  following  assumptions  are  made  in  order  to  simplify  the  radiation  heat  exchange  model 
and  to  make  its  calculation  tractable. 


iso-thermal  Each  gas  layer  and  each  wall  segment  is  assumed  to  be  at  a uniform 

temperature.  This  assumption  breaks  down  where  wall  segments  meet. 

equilibrium  The  wall  segments  and  gas  layers  are  assumed  to  be  in  a quasi-steady 

state.  The  wall  and  gas  layer  temperatures  are  assumed  to  change 
slowly  over  the  duration  of  the  time  step  of  the  associated  differential 
equation. 


fire  source  The  fire  is  assumed  to  radiate  uniformly  in  all  directions  giving  off 

a fraction,  x.  of  the  total  energy  release  rate  to  thermal  radiation. 
This  radiation  is  assumed  to  originate  from  a single  point.  Radiation 
feedback  to  the  fire  and  radiation  from  the  plume  is  not  modeled.  The 
plume  could  be  modeled  as  a collection  of  point  source  fires^,  however. 

®RAD4  allows  up  to  ten  point  source  fires 
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radiators 


opacity 


geometry 


The  radiation  emitted  from  a wall  surface,  a gas  and  a fire  is  assumed  to 
be  diffuse  and  gray.  In  other  words,  the  radiant  fluxes  emitted  by  these 
objects  are  independent  of  the  direction  and  the  wavelength.  They  can 
depend  on  temperature,  however.  Since  both  the  emittances  and  the 
temperatures  of  wall  segments  are  inputs  to  the  radiation  algorithms, 
it  is  assumed  that  the  emittances  are  consistent  with  the  corresponding 
wall  temperatures.  These  assumptions^  allow  us  to  infer  that  the  emit- 
tance,  e,  absorbance,  a,  and  reflectance,  p,  are  related  via  e = a = 1— p. 
A discussion  of  this  assumption  can  be  found  in  [12,  p.  589-590]. 

The  wall  surfaces  are  assumed  to  be  opaque.  When  radiation  encoun- 
ters a surface  it  is  either  reflected  or  absorbed.  It  is  not  transmitted 
through  the  surface.  Equation  (5)  would  have  to  be  modified  to  account 
for  the  loss  (or  gain)  of  energy  through  semi-transparent  surfaces. 

Rooms  or  compartments  are  assumed  to  be  rectangular  boxes.  Each 
wall  is  either  perpendicular  or  parallel  to  every  other  wall.  Radiation 
transfer  through  vent  openings,  doors,  etc  is  ignored. 


2.1.2  Derivation  of  the  Net  Radiation  Equations 

To  put  the  general  net  radiation  equations  in  perspective,  it  is  first  shown  how  to  compute 
the  radiation  exchange  between  N black  body  surfaces.  A general  enclosure  with  N walls 
containing  a transparent  gas  is  depicted  in  Figure  1.  Each  wall,  k,  has  a temperature  Tk 
and  an  area  Ak  ■ If  each  wall  segment  is  a black-body  then  reflection  terms  can  be  ignored 
when  computing  the  net  radiation  given  off  by  each  segment.  This  net  radiation  is  the 
difference  between  the  energy  given  off  and  the  energy  received  at  a wall  segment  as  in 


N 


AkAq"k  = O’ 


A,T}  - 

j = l 


where  is  the  fraction  of  energy  given  off  by  wall]  that  is  intercepted  by  segment  k.  The 
" notation  is  used  to  denote  flux,  (ie  flow  rate  of  mass  or  energy  per  unit  area).  The  above 
equation  can  be  simplified  using  the  symmetry  relation^,  AkFk-j  = AjFj^k,  to  obtain 


For  the  case  of  two  infinite  parallel  plates,  {N  = 2,  ^2-1  = 1))  and  F2-2  = 0 the  above 
equation  reduces  to  the  familiar  formula  for  radiation  exchange  between  two  bodies. 


Aq\  = ,7(7? -7?), 

= ,7(7?  - T?)  = -A?",  . 

■‘Diffusivity  implies  that  e;^  = a\  for  each  wavelength  A while  the  gray  gas/surface  assumption  implies 
that  is  constant  for  all  wavelengths 

^This  cind  other  configuration  factor  properties  will  be  discussed  in  section  2.1.5 
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For  the  more  general  case,  we  must  consider  wall  segments  that  are  not  black  bodies  and 
gas  layers  that  are  not  transparent.  For  this  case,  the  k’th  surface  emits  an  energy  flux  of 
crekT^  where  ct  < 1 and  a reflectance  term  consisting  of  radiative  fluxes  emitted  by  other 
surfaces.  The  radiation  exchange  at  the  k’th  surface  is  diagrammed  in  Figure  2.  For  each 
wall  segment  k from  1 to  AT  we  must  find  a heat  flux,  , such  that 


source  reflected  incoming  net 

+ \l = . 

Radiation  exchange  at  each  wall  segment  has  surface,  reflected,  incoming  and  net  radiation 
terms.  The  above  equation  then  represents  a system  of  linear  equations  that  must  be  solved 
for  Aq”  to  determine  the  net  fluxes  given  off  by  each  surface.  The  set  up  and  solution  of 
this  linear  system  is  the  bulk  of  the  work  required  to  implement  the  net  radiation  method 
of  Siegel  and  Howell. 

The  outgoing  radiation,  , consists  of  two  components,  gray-body  surface  radiation 
and  in-coming  radiation  from  other  surfaces  that  is  reflected.  The  incoming  radiation, 
is  composed  of  gray-body  surface  radiation  from  other  surfaces  in  the  enclosure,  radiating 
point-source  fires  and  emission  from  two  gas  layers.  The  radiation  leaving  the  k’th  surface 
is  given  by 


qr*  =Ak<T€kn  + il-ek)qi’^  . (1) 

The  radiation  arriving  at  the  k’th  surface  is 

N 

= (2) 
j=i 

where  is  a configuration  factor  and  is  defined  to  be  the  fraction  of  energy  leaving  a 

surface  j that  is  intercepted  by  surface  k.  The  term  r,_jt  is  a transmission  factor,  the  fraction 
of  energy  that  is  transmitted  unimpeded  through  the  gas  layer (s)  between  surfaces  j and  k. 
The  term,  ct,  represents  contributions  due  to  other  heat  source  such  as  heat  emitting  gas 
layers  and  point  source  fires.  The  computation  of  Fj^k  and  Tj_k  is  discussed  in  sections 
2.1.5  and  2.1.7  respectively.  The  computation  of  Ck  is  discussed  in  section  2.1.3. 

The  net  radiation,  AkAq'\,  is  the  difference  between  the  out-going  and  in-coming 
radiation  terms.  It  is  the  amount  of  energy  that  must  be  supplied  to  the  surface  k in  order 
to  maintain  a steady  state,  i.e.  constant  temperature  wall  segment.  This  is  the  net  radiation 
that  is  being  solved  for.  It  is  given  by 

AkAq\  = qr  - gr  (3) 

and  is  described  in  [7,  Chapters  8 and  17]. 

Equations  (l)-(3)  involve  equations  for  the  unknowns  gj.’^  and  Aq'\  for  k = 
1,  • • - , TV.  Siegel  and  Howell  [7]  find  the  net  radiant  flux  Aq" k for  each  of  N surfaces  in  an 
enclosure  by  solving  equations  (1),  (2)  and  (3)  for  Aq" k to  yield  the  net  radiation  flux.  The 
net  radiation  method  was  first  developed  by  Hottell  [8]. 

It  should  be  noted  that  the  net  radiation  method  (solving  for  Aq" f.)  is  superior  numeri- 
cally to  the  alternative  of  solving  equations  (l)-(3)  for  q'^^ . Though  equivalent  analytically, 
this  “new”  algorithm  is  inferior  numerically  to  the  net  radiation  method  (solving  for  Aq'\) 
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Figure  1:  N-Wall  Black  Body  Radiation  Exchange 
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Figure  2:  Energy  Distribution  at  the  k’th  Surface 


6 


because  of  the  catastrophic  cancellation  errors  introduced  by  the  subtraction  of  nearly  equal 
large  numbers  in  equation  (3).  The  same  comment  applies  when  solving  for  gj.". 

To  solve  by  the  net  radiation  method  for  Aq'\,  first  eliminate  gj.”  in  equations  (1)  and 
(2)  by  using  equation  (3)  to  obtain 


-out  _ A I ^^4  1 A 


qr  = 


■Ag' 


(4) 


N 


out 

<lk 


-AkAq\  = + . 

j=i 


The  term  g|“‘  can  be  eliminated  from  the  above  equation  by  using  (4)  to  obtain 


AktrT^  - Ak^^  = T 
i=i  ^ 


^ — —AjAqj  ) Fj.kTj-k  + Ck 

ft 


The  area  Ak  can  be  eliminated  after  rearranging  terms  and  using  the  symmetry  rela- 
tionship AkFk-j  = AjFj_k,  to  convert  the  above  equation  into  the  flux  equation 


Aq" 


N 


- E ^ , (5) 

J=1  ^ j=l 

Equation  (5),  listed  in  [7]  as  equation  (17-20)  is  the  net  radiation  equation.  The 
algorithms  used  in  this  report  use  equations  that  are  equivalent  analytically  to  the  net 
radiation  equations  but  different  (and  superior)  numerically.  These  modifications  and  their 
importance  are  discussed  in  section  3.1 
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2.1.3  Heat  Flux  Striking  a Wall  Segment 

In  general,  every  possible  path  between  two  wall  segments  should  be  considered  in  order 
to  compute  the  total  radiant  heat  transfer  between  these  segments.  This  is  not  practical 
in  a zone  fire  model  due  to  the  excessive  computational  costs.  The  approach  taken  here 
is  to  model  this  heat  transfer  using  just  one  path.  For  a typical  path  there  are  four  cases 
to  consider.  A path  from  wall  segment  j to  k can  start  in  either  the  upper  or  lower  layer 
and  may  finish  in  either  the  upper  or  lower  layer.  A fraction,  a = 1 — r,  of  the  energy 
encountering  a layer  is  absorbed.  The  rest,  r,  passes  through  unimpeded.  Table  1 gives 
formulas  for  the  heat  flux  striking  the  k’th  wall  segment  due  to  point  source  fires  and  heat 
emitting  gas  layers.  These  formulas  are  components  of  c^'  which  appear  in  equation  (23). 
Subsequent  sub-sections  discuss  how  to  compute  the  components  of  the  equations  in  Table 
1. 


Heat  Flux  Striking  a Wall  Segment  Due  to  a Point  Source  Fire  The  solid  angle 
of  a surface  is  the  fraction  of  the  total  radiant  energy  released  by  the  fire  that  strikes  the 
given  surface.  The  fire  is  assumed  to  radiate  uniformly  in  all  directions.  In  the  case  of  a 
point  source  a solid  angle  is  a configuration  factor.  If  the  gas  layers  are  transparent  then 
the  flux  striking  the  k’th  surface  due  to  the  f’th  fire  is 
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Table  1:  Radiative  Heat  Flux  Striking  the  k’th  Rectangular  Wall  Segment 


Path 

Fire 

Gcis  Layer 

9 f-k 

„L,gas 

9 j-k 

,,U,gas 

9 j-k 

in  upper 

-U 

^f-k  4TXfc 

0 

from  upper  to  lower 

Fk-j(Taf_kTl 

Fk-j(ray_kT^Tj^_,^ 

from  lower  to  upper 

_L  -U 

Fk-j<Taf_^TlTf‘_i^ 

Fk-i<raV_,Ti 

in  lower 

-L  xqVo'th'^j-k 
'f-k  4T.4fc 

Fk-j<Taf_f.Tl 

0 

_iljire  _ X^total‘^f->= 

’ ~ ■ 

The  total  energy  release  rate  of  the  fire  is  qlHh,  X is  the  fraction  of  this  energy  that 
contributes  to  radiation,  u;/_i/(47r^jfc)  is  the  fraction  of  the  radiant  energy  leaving  the  f’th 
fire  that  is  intercepted  by  the  k’th  wall  segment,  ie  a configuration  factor.  On  the  other  hand, 
if  the  gas  layers  are  not  transparent  then  there  are  four  cases  to  consider  when  calculating 
the  heat  transfer  from  a fire  to  a wall  segment.  The  fire  can  be  in  the  upper  or  lower  layer 
and  the  surface  can  be  in  the  upper  or  lower  layer.  Figure  3 shows  how  radiation  from  a 
fire  is  absorbed  by  each  layer  when  the  fire  is  in  the  lower  layer  and  the  surface  k is  in  the 
upper  layer.  The  other  three  cases  are  handled  similarly.  These  four  cases  are  summarized 
in  the  first  column  of  Table  1.  This  column  gives  formulas  for  the  flux  striking  a surface  k 
due  to  a point  source  fire. 

Heat  Flux  Striking  a Wall  Segment  Due  to  an  Emitting  Gas  Layer  The  energy 
emitted  by  the  i’th  layer  (i:=upper,  or  i=lower)  along  the  j-k’th  path  is 


9 3-k 


= ot)_^(TTi^ 


where  = 1 — The  emittance  of  the  gas  in  this  equation  is  the  same  as  the 

absorptance  due  to  the  gray  gas  assumption.  The  computation  of  a and  r is  discussed 
in  section  2.1.7.  Again  four  cases  must  be  considered  to  calculate  the  flux  striking  a wall 
segment.  The  last  column  of  Table  1 gives  formulas  for  radiation  striking  the  k’th  wall 
segment  due  to  gas  layer  heat  emissions  for  each  possible  path. 


2.1.4  Gas  Absorbance 

The  net  radiation  leaving  each  surface,  k,  is  k-  If  the  fire  is  ignored  then  the  total  energy 
added  to  the  gas  is 


N 
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Table  2:  Radiant  Heat  Absorbed  by  the  Upper  Layer 


Path  through  the 
Gas 

Due  to  Heat  Emitting 
Wall  Surface 

g°“fc  = AjFj-k 

Due  to  Gas  Layer 
Emission 

g")’lV=ct)_,<rn 

9]-V 

q"\'iVA,F,.k 

Due  to  Point  Source 
Fire 

9 f-k 

from  the  upper  to 
either  the  lower  or 
upper  layer 

-Out 

gj-k^^j-k 

„U,gaa 

9j-k 

from  the  lower  to 
the  upper  layer 

out  L U 

*lj-k^j-k°‘j-k 

L,gaa  U „U,gaa 

^J  — k ^]  — k ^]  — k 

from  the  lower  to 
the  lower  layer 

0 

0 

0 

If  the  gases  are  transparent  then  the  above  equation  must  sum  to  zero.  This  is  a good 
check  on  RAD2,  RAD4  and  RADIO.  The  energy  absorbed  by  the  gas  layers  may  be  due 
to  radiating  wall  segments,  emission  from  other  gas  layers  and  radiation  from  fires.  Tables 
2 and  3 summarize  the  formulas  used  to  compute  gas  layer  energy  gain/loss  due  to  these 
phenomena.  There  are  again  four  cases  to  consider,  since  an  arbitrary  path  may  start  in 
either  the  lower  or  the  upper  layer  and  end  in  the  lower  or  upper  layer.  Figure  4 illustrates 
the  heat  absorbed  by  the  gas  layers  due  to  surface  rectangle  emission  where  the  “from”  wall 
segment  is  in  the  upper  layer  and  the  “to”  wall  segment  is  in  the  lower  layer.  The  other 
three  cases  are  handled  similarly. 

2.1.5  Configuration  Factor  Properties 

A configuration  factor,  Fj^k,  is  the  fraction  of  radiant  energy  leaving  a surface  j that  is 
intercepted  by  a surface  k.  The  terms  used  to  define  a configuration  factor  are  illustrated  in 
Figure  5.  Configuration  factors,  also  called  view  or  shape  factors,  are  derived  in  numerous 
texts  [13,  Chapter  2],  [12,  Chapter  13],  [7,  Chapter  7].  This  sub-section  lists  some  configura- 
tion factor  properties  used  later  to  reduce  the  number  of  required,  direct  configuration  factor 
calculations  and  therefore  to  improve  the  efficiency  of  the  radiation  algorithm.  Other  con- 
figuration factors  are  calculated  indirectly  from  these  formulas  using  algebraic  relationships 
based  on  equations  (7)  - (11) 

The  configuration  factor  from  wall  segment  j to  wall  segment  k is  given  by 
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Figure  3:  Energy  Distribution  Due  to  a Fire 


leaving  surface  j q. 


out 


deposited  in  qout  ^ 
upper  layer  j J-k 


arriving  at  qOut  u 
interface  j-k 


deposited  in  qout^u 
lower  layer  j-k  j-k 


arriving  at  qOut  ^ u x !- 
rectangle  j j-k  i'*' 


Figure  4:  Energy  Distribution  Due  to  Wall  Segment  Emission 
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Table  3:  Radiant  Heat  Absorbed  by  the  Lower  Layer 


Path  through  the 
Gas 

Due  to  Heat  Emitting 
Wall  Surface 

qrJk=A,F,.,  {aTj- 

Due  to  Gas  Layer 
Emission 

q'r-V=<^)-k<^Tt 

9]-V 

/ti,gasA  p 

g j-fc 

Due  to  Point  Source 
Fire 

„///*re  _ 

? /-fc  - 

from  the  lower  to 
either  the  lower  or 
upper  layer 

L,gaa 

///ire  L 
q f-k^^f-k 

from  the  upper  to 
the  lower  layer 

q°'L\r^-kCtj-k 

_ „L,gaa 

^j  — k — k ^j  — k 

from  the  upper  to 
the  upper  layer 

0 

0 

0 

Figure  5:  Configuration  Setup 
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cos{6j)cos(0k) 


dAkdAj  . 


iA, 

The  following  properties  can  be  derived  from  this  definition  for  any  wall  segments 


(6) 


AjFj.k 

= AkFk-j 

(7) 

Fi-jQk 

= Fi-j  + F,_fc 

(8) 

i^jFiQj-k 

= AiFi^k  + -^jFj-k 

(9) 

where  i 0 j denotes  the  union  of  two  wall  surface  i and  j. 

If  the  collection  of  N wall  segments  form  an  enclosure  (he.  a room)  then  the  configuration 
factors  also  satisfy  equation  (10).  This  occurs  since  the  energy  leaving  the  j’th  surface  must 
strike  one  of  the  N surfaces  in  the  enclosure  (possibly  itself). 


N 

^Fj_k  = l j = (10) 

k=l 

If  four  wall  segments  1 • • - 4 are  configured  as  illustrated  in  Figure  6 then  the  following 
configuration  factor  identity  also  holds. 


= A2F2-3  (11) 

By  assuming  that  two  wall  segments  are  either  parallel  or  perpendicular  then  only  two 
classes  of  configuration  factors  need  to  be  calculated  directly  in  a zone  fire  model  radiation 
algorithm.  Equation  (6)  is  used  to  derive  configuration  factors  for  simple  geometries  that 
occur  in  zone  fire  modeling.  A catalog  of  configuration  factor  formulas  for  various  geometric 
configurations  is  given  in  [7,  Appendix  C].  Direct  configuration  factor  formulas  used  in 
RAD2,  RAD4  and  RADIO  follow. 

Consider  two  finite  rectangles  perpendicular  to  each  other  having  a common  edge  of  the 
same  length  as  illustrated  in  Figure  7.  The  dimensions  of  the  “from”  rectangle  is  / x it; 
and  the  dimension  of  the  “to”  rectangle  is  / x /i;  1 is  common  to  both  rectangles.  The 
configuration  factor,  <j)perp{h,l,w),  from  rectangle  1 to  2 is  then 


F\  — 2 — <^perp(h,  /,  l/t)  — 


1 


(lTTan-^47  + /fTan-i4  - V^^Tl^Tan-^ 

{ W H 


1 


+ 


1 / (1  + IT2)(1  + r ^2(1  + ^ 

4 \ + + [(1  + IT2)(//2  + ^2) 


+ + W^) 


(1  + + IT2) 


(12) 


where  H = h/ 1 and  W = w/ 1 . This  formula  is  found  in  [7,  p.  825]. 

We  also  need  to  compute  the  configuration  factor  between  two  identical,  parallel,  directly 
opposed  rectangles.  This  situation  is  setup  in  Figure  7.  The  distance  between  the  two 
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=^2"2-3 


Figure  6:  Configuration  Factor  Symmetries 


a 


Figure  7:  Configuration  Factor  Setup  For  Pairs  of  Parallel  and  Perpendicular  Rectangles 
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rectangles  is  c.  The  length  and  width  of  these  rectangles  is  a and  6.  The  configuration 
factor,  <t>par{‘^ib,c)  between  these  two  rectangles  is  then 


Fi  — 2 — 4>pario,  b,  c)  — 


+ y2 


+ 


y\/l  + X2Tan-i 


1 + X2 


(13) 


where  X = a/c  and  Y = bfc.  This  formula  is  symmetric  in  X and  Y.  It  can  be  found  in  [7, 
p.  824], 

For  a room  with  N wall  segments,  N x N = configuration  factors  must  be  calculated. 
Equations  (12)  and  (13)  are  expensive  to  compute  due  to  the  complicated  expressions  in- 
volving log  and  Tan~^  functions.  This  portion  of  the  work  is  reduced  in  RAD4  by  noting 
that  only  2 configuration  factor  calculations  involving  equation  (12)  are  required  rather  than 
4 X 4 = 16.  The  other  14  configuration  factors  are  obtained  using  algebraic  relationships. 
These  algebraic  formulas  are  given  in  section  2.2.2.  For  the  RADIO  case  only  eight  con- 
figuration factor  calculations  need  be  calculated  using  equations  (12)  and  (13).  Again,  the 
other  92  can  be  obtained  using  algebraic  formulas.  These  formulas  are  detailed  in  section 
2.2.3. 


2.1.6  Solid  Angles 

A solid  angle  as  illustrated  in  figure  8 is  the  area  intercepted  on  a unit  sphere  by  a conical 
angle  originating  at  the  sphere  center.  A solid  angle  is  used  to  determine  the  fraction  of  a 
radiating  point  source  fire  that  strikes  a surface.  This  formula  was  obtained  from  Yamada’s 
formula  in  [10]  by  observing  that 


Sin 


-1 


y/x^  + y^ 


-I-  Sin 


-1 


y/x^  + y^ 


The  solid  angle  of  a rectangle  with  sides  of  length  x and  y that  lies  in  a plane  a distance 
r from  the  center  of  a sphere  is  given  by 


= AJ  Sin 


-1 


\/y^  + 


-f-  Sin 


-1 


y/  x^  + 


(14) 


where 


The  solid  angle  uj(x,  y)  is  symmetric  in  x and  y.  Solid  angles  are  also  additive,  so  that 
the  solid  angle  of  an  arbitrary  rectangle  can  be  computed  using  (14)  and 


^(a;i,X2,yi,y2)  = Sgn(a:2y2)w(|x2|,  |y2|)  - Sgn(a;2yi)t^(k2|,  |yi|) - 

sgn(xiy2)w(Ixi|,  |y2|)  -|-  sgn(xiyi)u;(lxi|,  |yi|) 
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0^  = solid  angle  for  rectangle  i 


CQ4  = COi+2+344  - COi+2  - Wi+3  + COi 


Figure  9:  Solid  Angles  For  Arbitrary  Rectangles 
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sgn(x)  = < 


1 if  X > 0 
-1  ifx<0. 


Figure  9 illustrates  how  the  solid  angle  of  rectangle  4 can  be  computed  in  terms  of  solid 
angles  of  rectangles  1©2©30  4,  1©3,  1©2  and  1. 


2.1.7  Transmission  Factors 

A transmission  factor,  r,  is  the  fraction  of  energy  passing  through  a gas  unimpeded.  The 
transmittance  of  a gas  depends  on  the  absorption  coefficient  of  the  gas  and  the  length  of 
the  path  through  the  gas.  A simple  relationship  for  r can  be  determined  by  assuming  that 
the  absorptance  of  the  gas  (a  local  phenomena)  is  uniform  throughout  the  gas  layer.  This 
factor,  a decaying  exponential,  is  given  by 

t{L)  = 

where  A is  the  absorptance  of  the  gas  per  unit  length  and  L is  the  path  length.  This  formula 
is  known  as  Beer’s  law.  The  gas  absorptance  is  not  calculated  by  the  radiation  exchange 
algorithms  presented  in  this  report.  Modak  in  [14]  gives  an  algorithm  for  calculating  gas 
absorptance  from  such  information  as  soot  concentration,  partial  pressures  of  CO,  C02  etc. 
The  emittance  of  the  gas  is  the  same  as  its  absorption  due  to  the  gray  gas  assumption.  The 
transmission  factor,  r in  the  above  equation  is  defined  for  one  specific  path  through  a gas.  We 
are,  however  considering  radiation  exchange  between  a pair  of  finite  area  rectangles  where 
many  paths  of  different  lengths  occur.  Siegel  and  Howell  define  an  average  transmission 
factor  [7,  p.  603]  considering  all  possible  paths  between  two  surfaces  through  the  gas.  This 
form  of  T is  defined  to  be 


-IL 


A,Ak 


t{L)  cos{9ic)  cos{9j) 


dAj  dAk/{AjFj^k) 


The  numerator  of  the  fraction  in  in  the  above  equation  is  the  same  as  a configuration 
factor  if  r(L)  = 1.  It  is  not  practical  to  compute  a transmission  factor  using  this  integral. 
We  can  estimate  this  integral  by  finding  a characteristic  path  with  length  L such  that 


Tj-k  = e 


-AL 


For  the  ten-wall  model,  this  path  is  taken  to  be  between  the  centers  of  two  rectangles. 
This  length  is  an  underestimate  of  L.  This  approximation  breaks  down  when  the  two  wall 
segments  are  close  together  or  one  of  the  wall  segments  is  a complex  shape  (such  as  the 
union  of  four  upper  walls).  The  four  and  two-wall  model  estimates  L based  upon  an  average 
distance  between  the  rectangles  that  make  up  the  wall  segments. 

For  a given  path  between  surface  j and  surface  k we  need  to  calculate  the  path  length 
through  the  upper  layer,  Lu  , and  the  path  length  through  the  lower  layer  L^.  Transmission 
factors  for  the  upper  and  lower  layers  are  then  defined  to  be 


rU  _ .-L^.^Au 


‘j-k 


= e 


^i-k  — ^ ^ 
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The  fraction  of  energy  that  passes  through  both  layers  is  then 

— k^j  — k • 

The  amount  of  energy  absorbed  by  a layer  is  just  the  amount  that  doesn’t  pass  through  a 
layer  or 


0‘f-k 


1 - = 1 - , 

1 — Tjlfc  = 1 — . 


2.2  Two,  Four,  Ten  Wall  Segment  Radiation  Exchange 

The  formulas  in  section  2.1  for  computing  radiation  exchange  were  specified  in  terms  of 
general  wall  segments.  This  section  discusses  the  radiation  exchange  computation  in  terms 
of  a two- wall,  four- wall  and  ten-wall  model. 


2.2.1  Two- Wall  Configuration  Factors 

The  two-wall  model  combines  the  ceiling  and  four  upper  walls  into  one  wall  segment  and 
the  four  lower  walls  and  the  floor  into  the  second  wall  segment.  The  configuration  factors 
for  these  two  surfaces  are  derived  by  Quintiere  in  [15,  Appendix]  and  are 


Fi_i 

Fi_2 


F2-1 


F2-2 


’ 

A2  ’ 


where  Ai,  Ad  and  A2  are  the  areas  of  the  extended  ceiling,  layer  interface  and  extended 
floor  respectively.  These  configuration  factors  are  used  in  the  original  two-wall  radiation 
model  in  CFAST  [5]  and  in  BRI  [3,  4]. 

The  two- wall  model,  RAD2,  interacts  with  a four-wall  heat  conduction  model  in  CFAST. 
The  ceiling  and  upper  wall  temperatures  may  be  different,  so  the  question  of  how  to  represent 
the  extended  ceiling  temperature  arises.  RAD2  chooses  an  extended  ceiling  temperature  that 
results  in  the  same  energy  contribution  to  the  enclosure  that  a four-wall  radiation  algorithm 
would  predict.  The  energy  added  to  the  room  due  to  the  ceiling  and  upper  wall  temperatures 
of  Tia  and  Tu  is 


{Aiaf-loTia  + AiJfiftTjj)  . 

We  want  to  choose  an  effective  or  average  temperature,  Tavg,  and  emittance,  €avg  for  the 
extended  ceiling  that  matches  this  energy  contribution. 

energy  from  extended  ceiling  energy  from  ceiling  energy  from  upper  wall 

/ . / ^ 

{Ala  + An)  (avgT^vg  = AiaClaT^a  + (1^) 
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An  average  emittance  for  the  two  wall  segments  is  computed  using  an  average  of  €ia  and 
eii  weighted  by  wall  segment  areas  or 


^avg 


A\a(.\a  + 

Ala  + Alb 


Equation  (15)  can  now  be  solved  for  Tavg  using  this  value  of  Cavg  to  obtain 


• avg 


AlgClgT^^  + AlbClbT^l, 
-dla^la  + Aib€ib 


A similar  procedure  can  be  used  to  compute  an  effective  temperature  and  emittance  for  the 
extended  floor. 


2.2.2  Four- Wall  Configuration  Factors 

The  configuration  factors  for  four-wall  radiation  exchange  are  derived  similarly  to  Quintiere’s 
derivation  for  two  walls  in  [15,  Appendix].  The  setup  for  the  following  derivation  is  given 
in  Figure  11.  We  wish  to  determine  the  configuration  factors 


F,_j  for  2,  j = 1,  - . . , 4. 

The  16  configuration  factors  can  be  determined  in  terms  of  Fi_4,  Fi^j  and  F^_d-  Fi_4 
does  not  change  during  a simulation  since  its  value  depends  only  on  the  height  of  the  room 
and  the  area  of  the  floor.  Therefore,  Fi_4  only  needs  to  be  computed  once.  Configuration 
factors,  Fi^d  and  F^^d  depend  on  the  layer  interface  height  so  need  to  be  calculated  each 
time  the  radiation  exchange  is  to  be  calculated.  Configuration  factors  Fi_4,  Fi^d  and  F<i^d 
are  determined  using  equation  (13).  Since  Ai  = A^  it  follows  that  F4_i  = Fi_4.  The  other 
14  configuration  factors  can  be  calculated  using  simple  algebraic  formulas. 

Since  the  floor  and  the  ceiling  is  assumed  to  be  a flat  rectangular  surface  it  follows  that 

Fi_i  = F4_4  = 0. 


Using  the  fact  that  configuration  factors  in  an  enclosure  sum  to  1 and  that  due  to  symmetry 
F2-1  = F2-d,  it  follows  that 


Fi_2  -b  Fi_(i 


1, 


F2-1  + F2-2  + F2-d  = 2F2_i  -f-  F2-2  — 2-— Fi_2  -t-  F2-2  — 1 

A2 


Similarly, 


solved  for 

Fi_ 

.2  and 

Fi_2  = 

1 - 

Fi.d 

F2-I  = 

A.2 

Fi-2 

F2-2  = 

1 - 

2F2-1 

F4-3  = 

1 - 

Fi-d 

II 

1 

A4 

-^3 

1 

U 

II 

CO 

1 

1 - 

2F3_4 

(16) 

(17) 
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Surface  1,  The  extended  ceiling 
surface  d,  the  layer  interface 


Surface  2,  The  extended  floor 


Figure  10:  Extended  Floor,  Ceiling  Configuration  Factor  Setup 


Surface  1,  The  ceiling 

Surface  2,  The 
upper  walls 

Surface  d,  The  layer  interface 


Surface  3,  The 
lower  walls 

Surface  4,  The  floor 


Figure  11:  Four  Wall  Configuration  Factor  Setup 
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Using  the  above  configuration  factors  and  equation  (10)  it  follows  that 


^1-3 

^3-1 

F3-2 

■P’2-3 

-P’2-4 

-P’4-2 


1 — Fi_4  — Fi_2 


1-F3_i 


F3-3  — F3_4 


1 — F2_i  — F2-2  — F2-3 


2.2.3  Ten-Wall  Configuration  Factors 

To  handle  the  more  general  radiation  exchange  case,  a room  is  split  into  ten  surfaces  as 
illustrated  in  Figure  12.  These  surfaces  are  the  ceiling,  four  upper  walls,  four  lower  walls 
and  the  floor.  The  radiation  exchange  is  computed  between  these  ten  surfaces  and  the 
intervening  gas  layer(s).  In  general,  100  configuration  factors,  Fj_jt  and  100  transmission 
factors  rj_jfe  need  to  be  determined  each  time  this  algorithm  is  invoked.  Although  there 
are  100  configuration  factors  for  this  room,  only  eight  have  to  be  calculated  directly  using 
equations  (12)  and  (13).  The  other  92  can  be  computed  in  terms  of  simple  algebraic  rela- 
tionships using  the  properties  outlined  in  equations  (7)  to  (11).  This  reduction  in  required 
configuration  factor  calculations  is  due  to  the  fact  that  the  rectangle  pairs  2 and  4,  3 and 
5,  6 and  8 and  7 and  9 each  have  equal  areas. 

The  following  15  configuration  factors  are  computed  only  once  during  a fire  simulation, 
assuming  that  the  room  does  not  change  size.  These  factors  are  between  the  six  surfaces  that 
form  the  enclosure.  There  are  a total  of  36  configuration  factors  between  these  six  surfaces. 
The  six  factors,  Fk-k,  are  each  zero.  Fifteen  factors  are  given  below.  The  other  15  can  be 
derived  using  the  symmetry  relation  AjFj^k  = AkFk-j.  The  first  six  configuration  factors 
are  computed  using  formulas  (12),  (13).  The  other  nine  configuration  factors  are  derived 
from  these  six  using  the  identities  (7)  through  (11).  These  formulas  are  implemented  in  the 
FORTRAN  subroutine  RMFIG 


Fi_2©6 

-P’1-307 

Fi_io 

^206-307 

F206-,408 

^307-509 


^perp  (Az,  Ax,  Ay) 
*/*perp(Az,  Ay,  Ax) 
4^par  (Ax,  Ay,  Az) 
<^perp(Ay,  Az,  Ax) 
(f^par  (Ax,  Az,  Ay) 
4^par  (Az,  Ay,  Ax) 


where  <i)perp  is  the  configuration  factor  for  two  perpendicular  rectangles  sharing  a common 
edge  given  by  equation  (12)  and  <f)par  is  the  configuration  factor  for  two  identical  parallel 
rectangles  given  by  equation  (13). 


Fi_408  — Fi_206 
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2 = upper  front 


8 = lower  back 


4 = upper  back 


1 = top 


10  = bottom 


6 = lower  front 


Figure  12:  Ten  Wall  Configuration  Factor  Setup 
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The  following  configuration  factors  are  computed  each  time  the  radiation  exchange  is 
computed  since  these  terms  depend  on  the  layer  height  which  may  vary.  The  first  eight 
configuration  factors  are  computed  using  formulas  (12),  (13).  The  other  92  configuration 
factors  are  derived  from  these  eight  using  the  identities  (7)  to  (11).  The  formulas  are  listed 
in  the  order  that  they  are  computed  so  that  once  a configuration  factor  appears  on  the  left 
hand  side  of  an  equation  (such  as  Fi^e)  it  may  be  used  later  on  the  right  hand  side.  These 
formulas  are  implemented  in  the  FORTRAN  subroutine  LAYFIG. 
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Second  Row  and  Column 


i^2-2  = 0 
F2-5  = F2-3 

F2-6  = 0 

The  formula  for  ^2-7  will  be  derived.  Derivations  for  other  factors  follow  similarly.  Equa- 
tions (8)  and  (9)  can  be  used  to  simplify  vl2e6-p2®6-3®7  to 

'42®6-^2®6-3®7  = j42E2-3®7  + -<46-f’6-3®7 

= A2F2-3  ■^2F2-7 AeFe^s  + AeFQ_7  . (18) 

The  factors  F2-7  and  Fq^s  are  related  via  A2F2-7  = AqFq^z  since  the  surfaces  2,  3,  6 and 
7 are  configured  as  illustrated  in  figure  6.  Equation  (18)  can  then  be  simplified  to 

■42®6.^2®6-3®7  = A2F2-3  + ^A2F2-7  + AeFe^7 

Solving  for  F2-7  yields  the  desired  formula  as  listed  below.  Note  that  F2Q3-3Q7,  F2-3  and 
F6~7  were  computed  previously. 
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Fourth  Row  and  Column 


F4_4  = 0 

F4-5  = F2-3 

F4_6  = F2-8 

F4-7  = F2-7 

F4_8  = 0 

F4_9  = F2-9 

F4-IO  = -^2-10 
Fj-4  = ^F4-jj  = 5,...,10 

Fifth  Row  and  Column 

i^5-5  = 0 

F5-6  = ^3-6 
F5-7  = ^3-9 
F5-8  = ^3-8 
F5-9  = 0 

-fs-lO  = -^3-10 

Fj-^  = 4^F5_,i  = 6,...,10 

Sixth  Row  and  Column 
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Seventh  Row  and  Column 
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Ninth  Row  and  Column 
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Tenth  Row  and  Column 
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■P’10-10  = 0 

3 Solving  the  Net  Radiation  Equations 

3.1  Solving  The  Net  Radiation  Equations  Efficiently 

The  net  radiation  equation  (5)  is  not  diagonally  dominant.  Iterative  methods  then  should 
not  be  used  to  solve  this  equation  unless  it  is  suitably  transformed.  This  can  be  done  by 
substituting,  Aq'\  = CfcAqJ,'  into  equation  (5)  to  obtain 

N N 

i=l  i=l  ^ 

There  are  two  reasons  for  solving  equation  (19)  instead  of  (5).  First,  since  Cfc  does  not  occur 
in  the  denominator,  radiation  exchange  can  be  calculated  when  a wall  segment  emittance 
is  zero.  Second  and  more  importantly,  the  matrix  corresponding  to  the  linear  system  of 
equations  in  (19)  is  diagonally  dominant.  When  the  number  of  wall  segments  is  large  and 
the  wall  segments  have  emittances  close  to  one  which  often  occurs  in  typical  fire  scenarios, 
the  time  required  to  solve  this  modified  linear  system  can  be  significantly  reduced  due  to 
this  diagonal  dominance  by  using  iterative  methods. 

To  see  this,  re-write  equation  (5)  into  matrix  form  to  obtain 

.4Ag"  = BE-c  (20) 


where  the  components  of  the  N x N matrices  A and  B are 


Pfc.i  771  _ 1 ~ 

O/tJ  = tk-jTj-k 


bk,j  = 6k,j  — Fk-jTj^k 
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(21) 

(22) 


and  the  k’th  component  of  the  vectors  d'  and  E are 


N 

Nf.re 

Ck 

Ak 

_ („iiU,ga$  ,tL,gaa\ 

- i-k  +9  i-k  J 

+ E d/:i , 

(23) 

J = 1 

/=i 

Ek 

II 

(24) 

and  j is  the  dirac  delta  function,  -Ffc-;  is  the  configuration  factor  from  the  k’th  to  the  j’th 
wall  segment,  Cj  is  the  emittance  of  the  j’th  wall  segment,  is  a fraction  ranging  from 
0 to  1 indicating  the  amount  of  radiation  that  is  transmitted  through  a gas.  Also, 
and  radiation  contributions  due  to  the  gas  layers  and  q"j'^l  are  radiation  terms 

due  to  the  f’th  fire.  The  matrix  A can  be  transformed  into  a diagonally  dominant  matrix 
using  the  following  scaling  matrix, 

/ei  0 0 \ 

^=  0 •••  0 
Vo  0 

where  €k  is  the  emittance  of  the  k’th  wall  segment.  Define  the  scaled  matrix  A by  post- 
multiplying  A by  D and  pre-multiplying  Aq"  by  D~^  to  obtain 

A = AD 
Aq"  = D-^Aq". 

Equation  (20)  then  reduces  to 

AAq"  = AAq"  = BE- c. 

Once  the  solution  Aq"  is  found  we  may  recover  the  solution,  Aq" , to  the  original  problem 
by  using  Aq"  = DAq".  The  matrix  A is  diagonally  dominant  which  is  now  shown. 

Using  the  definition  of  Okj  in  equation  (21)  the  kj’th  element  of  A is 

dfc,j  — ^k,j^j  — ^k,j  — Tfc_j7j_fc(l  fj)-  (25) 

A matrix  is  diagonally  dominant  if  for  each  row  the  absolute  value  of  the  diagonal  element 
is  greater  than  the  sum  of  the  absolute  values  of  the  off  diagonal  elements  or  equivalently 

N 

iSfc.fcl  > ^ (26) 

j = 1 
ji^k 

Substituting  (25)  into  (26)  we  get  the  following  requirement  for  A to  be  diagonally  dominant. 

N 

^ — Ek-k{^  — £k)Tk-k  > ^ Fk-j  (I  — £j)  Tj^k 

j = 1 
j 
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or  equivalently 


N 

j=i 

The  matrix  A is  then  diagonally  dominant  since  I > {1  — ej)Tj_k  and 


N N 

1=  ^k-j>  XI  Fk-jil  - £j)Tj_k  . 

j = 1 j = I 


Iterative  techniques®  are  guaranteed  to  converge  for  diagonally  dominant  matrices  [16, 
p.  542].  They  also  can  be  much  more  efficient  than  direct  methods  such  as  Gaussian 
elimination.  The  convergence  speed  depends  on  how  small  the  right  hand  side  of  the  above 
inequality  is  compared  to  1.  Physically,  if  the  surfaces  being  modeled  are  approximate  black 
bodies  (e  close  1)  or  the  gas  layers  are  thick  (r  close  to  0)  then  iterative  techniques  for 
solving  the  net  radiation  equations  will  converge  rapidly.  Typical  emittances  for  materials 
used  in  fire  simulations  range  from  e = .85  to  .95.  For  the  limiting  case  when  the  wall 
materials  are  black  bodies  then  the  matrix  .A  is  a diagonal  matrix  and  iterative  methods 
will  converge  in  one  iteration. 

The  advantage  of  using  an  iterative  method  over  a direct  method  for  computing  radia- 
tion exchange  between  approximate  black  bodies  increases  as  the  number  of  wall  segments 
increases.  The  cost  of  solving  the  linear  system  directly  is  proportional  while  the  cost 
of  using  iterative  techniques  is  proportional  to  kN^  where  k is  the  number  of  iterations  and 
N is  the  number  of  wall  segments.  Using  Gauss-Seidel  iterative  methods,  it  has  been  found 
that  convergence  is  achieved  after  two  to  three  iterations  for  emittances  around  .9.  The 
break-even  point  between  iterative  and  direct  methods  for  matrices  of  size  10  is  about  6 
or  7.  The  linear  system  for  RAD2  and  RAD4  is  2x2  and  4x4  respectively.  Iterative  meth- 
ods are  not  faster  for  problems  this  small.  RADIO  and  problems  with  more  wall  segments 
can  use  iterative  methods  to  decrease  the  time  required  to  solve  the  linear  system  without 
sacrificing  accuracy. 

The  radiation  exchange  equations  can  be  solved  analytically  using  Cramer’s  rule  for  the 
two  wall  segment  case.  This  is  how  the  radiation  exchange  equations  were  derived  in  [3] 
and  [5].  Cramer’s  rule  is  not  a good  numerical  technique  to  use  for  the  solution  of  linear 
systems  (even  for  2x2  systems)  due  to  cancellation  error  that  can  be  introduced  when 
solving  equations  that  are  ill-conditioned. 

3.2  Algorithm  for  Calculating  Four  Wall  Radiation  Exchange 

The  strategy  for  computing  the  radiation  exchange  between  four  wall  segments  is  outlined 
below.  RAD4  performs  these  steps  directly  or  call  subroutines  that  perform  them.  RAD2 
and  RADIO  follows  the  same  logic. 

Input 

Temperatures  Ceiling,  Upper  Wall,  Lower  Wall,  Floor 
®Gauss-Seidel  or  for  example. 
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Emissivities  Ceiling,  Wall,  Floor 

Absorptivities  Upper,  Lower  Layer 

Fire  Size,  Location,  number 

Room  room  number,  dimensions,  layer  height 

Output 

Flux  ceiling,  upper  wall,  lower  wall,  floor 
Energy  Absorption  Rate  upper  layer,  lower  layer 

Steps  1.  Calculate  configuration  factors,  solid  angles  as  described  in  section  2.2.2  and 

2.1.6. 

2.  Determine  the  effective  length  between  each  pair  of  wall  segments.  From  these 
lengths  and  inputted  layer  absorptivities  calculate  transmission  factors  for  surface 
j to  surface  k 

3.  Calculate  transmission  factors  and  gas  layer  absorptions  for  each  fire  f to  surface 
k as  outlined  in  section  2.1.7. 

4.  Calculate  the  energy  absorbed  by  each  gas  layer  due  to  upper/lower  gas  layer 
emission  and  due  to  the  fire(s)  following  tables  2 and  3. 

5.  Set  up  the  linear  algebra 

(a)  Define  vector  E using  equation  (24) 

(b)  Define  matrix  A using  equation  (25) 

(c)  Define  matrix  B using  equation  (22) 

(d)  Define  vector  c,  using  equation  (23)  and  Table  1. 

6.  Solve  the  linear  system 

iA?"  = BE- c (27) 

for  Aq"  net  radiation  leaving  each  surface  where  Aq" f.  = Aq'^Ck  • If  the  emittances 
are  sufficiently  close  to  1 then  use  iteration  to  solve  equation  (27)  otherwise  use 
Gaussian  elimination.  The  reason  for  this  is  due  to  the  diagonal  dominance  of  A 
as  explained  in  the  previous  section. 

7.  Calculate  the  energy  absorbed  by  the  upper  and  lower  gas  layers  due  to  the  total 
energy  , leaving  each  rectangle  k. 


4 Computational  Results 

4.1  Checks 

Several  simple  checks  can  be  made  to  verify  a portion  of  the  radiation  calculation.  First,  no 
heat  transfer  occurs  when  all  wall  segments  and  both  gas  layers  are  at  the  same  temperature. 
Therefore,  the  net  radiation  flux,  Aq" ^ given  off  by  each  surface  and  the  energy  absorbed 
by  the  gas  should  be  zero  under  these  uniform  temperature  conditions.  Second,  when  there 
is  no  fire,  the  net  energy  absorbed  by  the  gas  must  be  the  same  as  the  net  energy  given  off 
by  the  wall  segments  or  equivalently 
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Table  4:  Four  Wall  Radiation  Algorithm  Timings 


Case 

Total  Time  (s) 

Linear  Solve  Time  (s) 

45  configuration  factors 
and  direct  linear  solve 

.2 

.06 

Eight  configuration  fac- 
tors and  iterative  linear 
solve 

.06 

.01 

Two  configuration  fac- 
tors and  direct  linear 
solve 

.012 

.001 

N 

qiower  + qupper  = energy  absorbed  by  interior  gases  = ^ Ak^q'\ 

k=l 

When  the  layers  are  transparent  then  the  above  equation  sums  to  zero  even  though  the 
individual  wall  fluxes  Aq'\  will  in  general  be  non-zero.  The  gas  absorbance  terms,  qiower 
and  qupper,  are  computed  by  RAD2,  RAD4  and  RADIO.  These  values  can  be  summed  to 
verify  that  above  equation  is  satisfied.  These  checks  were  found  to  be  useful  since  they  all 
failed  to  hold  during  some  point  in  the  development  process! 

4.2  Timings 

Configuration  calculations  are  one  of  the  major  bottle  necks  in  the  radiation  exchange 
calculation.  Techniques  to  reduce  the  number  of  these  calculations  will  improve  the  algo- 
rithms efficiency.  The  original  version  of  RAD4  was  based  on  RADIO,  a ten  wall  segment 
model.  It  computed  45  configuration  factors  directly.  Subsequent  versions  of  RAD4  com- 
puted eight  and  two  configuration  factors  directly.  Table  4 summarizes  the  time  required  by 
these  three  different  versions  of  RAD4.  This  table  shows  that  the  first  version  of  RAD4  used 
approximately  70%  of  the  time  to  setting  up  the  linear  system  and  30%  solving  it.  Reducing 
the  setup  overhead  by  computing  fewer  configurations  factors  reduced  the  computation  time 
required  by  a factor  of  17. 

Why  quibble  over  the  timings  of  a subroutine  that  only  took  .2  seconds  to  execute  to 
begin  with?  The  relative  impact  of  RAD4  on  CFASTs  performance  can  be  measured  by 
comparing  the  time  required  to  execute  DSOURC^  with  and  without  RAD4.  For  a six  room 
CFAST  test  case  DSOURC  took  about  .06  seconds®.  A routine  that  takes  .2  seconds  per 
room  used  in  each  room  will  result  in  a 21-fold^  increase  in  computer  time.  Even  the  fastest 
version  of  RAD4  will  cause  an  increase  of  execution  time  of  2.25  if  it  is  used  in  each  room. 

^DSOURC  is  the  subroutine  CFAST  uses  to  calculate  the  right  hand  side  of  the  modeling  differential 
equations.  Most  of  the  work  is  performed  by  this  routine  or  routines  that  DSOURC  CcJls. 

*A11  times  are  measimed  on  a Compciq  386/20  Deskpro.  This  computer  has  a 20mhz  clock  and  uses  a 
floating  point  accelerator  (math  co-processor).  The  actual  times  will  be  different  on  different  computers. 
But  the  relative  times  and  hence  the  conclusions  should  be  the  same. 

®DSOURC  time  with  RAD4/DSOURC  time  withoutRAD4  = (.2  • 6 -|-  .06)/. 06  a 21 
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4.3  Comparisons  of  RAD2  with  RAD4 

This  subsection  compares  the  predictions  of  a two  wall  radiation  exchange  model,  RAD2, 
with  a four-wall  model,  RAD4.  One  of  the  assumptions  made  in  section  2.1.1  about  N 
wall  segment  radiation  models  is  that  the  temperature  distribution  of  each  wall  segment  is 
approximately  uniform.  The  zone  fire  model  CFAST  models  the  temperature  of  four  wall 
segments  independently.  Therefore,  a two  wall  model  for  radiation  exchange  can  break  down 
when  the  temperatures  of  the  ceiling  and  upper  walls  differ  significantly.  This  could  happen 
in  CFAST,  for  example,  when  different  wall  materials  are  used  to  model  the  ceiling,  walls 
and  floor.  To  demonstrate  this  consider  the  following  example. 

To  simplify  the  comparison  between  the  two  and  four  wall  segment  models,  assume 
that  the  wall  segments  are  black  bodies  (the  emissivities  of  all  wall  segments  are  one) 
and  the  gas  layers  are  transparent  (the  gas  absorptivities  are  zero)  . This  is  legitimate 
since  for  this  example  we  are  only  interested  in  comparing  how  a two  wall  and  a four  wall 
radiation  algorithm  transfer  heat  to  wall  segments.  Let  the  room  dimensions  be  4 x 4 x 4 
[m],  the  temperature  of  the  floor  and  the  lower  and  upper  walls  be  300  [K].  Let  the  ceiling 
temperature  vary  from  300  [K]  to  600  [Kj. 

Figure  13  shows  a plot  of  the  heat  flux  striking  the  ceiling  and  upper  wall  as  a function 
of  the  ceiling  temperature.  The  two  wall  model  predicts  that  the  extended  ceiling  (a  surface 
formed  by  combining  the  ceiling  and  upper  wall  into  one  wall  segment)  cools,  while  the  four 
wall  model  predicts  that  the  ceiling  cools  and  the  upper  wall  warms.  The  four-wall  model 
moderates  temperature  differences  that  may  exist  between  the  ceiling  and  upper  wall  (or 
floor  and  lower  wall)  by  allowing  heat  transfer  to  occur  between  the  ceiling  and  upper  wall. 
The  two  wall  model  is  unable  to  predict  heat  transfer  between  the  ceiling  and  the  upper 
wall  since  it  models  them  both  as  one  wall  segment. 

A four-wall  algorithm  will  also  break  down  when  the  uniform  temperature  assumption  is 
broken.  This  could  occur  when  a fire  is  located  nearer  to  one  side  of  a room  than  another. 

5 Conclusions 

This  report  documented  algorithms  for  computing  radiative  heat  exchange  for  three  special 
cases,  a two-wall,  four-wall  and  ten-wall  model.  The  theoretical  basis  for  a general  N 
wall  model  is  well  documented  in  the  literature  [7,  8,  9].  But  an  implementation  of  an  N 
wall  model  is  not  yet  practical  for  a zone  fire  model  due  to  the  high  computational  costs 
compared  to  other  components  in  a zone  fire  model.  One  step  was  taken  towards  making  N 
wall  models  practical.  For  wall  surfaces  that  are  approximate  black  bodies  (e  > .85)  it  was 
shown  that  the  linear  system  of  equations  involving  the  unknown  net  radiative  flux  could 
be  solved  iteratively,  reducing  an  o{N^)  to  an  o(N^)  problem.  For  specific  cases,  (four-wall, 
ten-wall)  it  was  shown  how  to  set  up  this  linear  system  efficiently  by  avoiding  unnecessary 
configuration  factor  calculations.  For  the  general  N wall  problem  it  is  not  enough  to  solve 
the  linear  system  efficiently.  For  example,  in  the  ten-wall  case  the  linear  solve  time  is  only 
20%  of  the  total  solution  time.  Methods  need  to  be  found  to  calculate  configuration  factors 
more  efficiently,  perhaps  at  some  cost  in  accuracy. 
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Figure  13:  Ceiling,  Upper  Wall  and  Extended  Ceiling  Flux  vs.  Ceiling  Temperature 
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Nomenclature 


A area  [m^ 

A absorption  coefficient  [m“^] 

A coefficient  matrix  for  the  net  radiation  equation. 

A coefficient  matrix  for  diagonally  dominant  version  of  the  net  radiation 

equation 

c vector  of  source  terms  used  in  the  net  radiation  equation  (5)  to  repre- 

sent energy  contributions  to  wall  segments  due  to  gas  emitting  layers 
and  point  source  fires  [W] 

B matrix  used  on  the  right  hand  side  of  the  net  radiation  equation 

E emmissive  power,  a vector  whose  Ar’th  component  is  crT^ 

D diagonal  scaling  matrix  used  to  convert  A to  the  diagonally  dominant 

version  A 


Fj-k 

N 

^fire 

? 

Aq" 

T 


6 


€ 


X 


geometric  configuration  factor,  also  called  a view  factor.  The  fraction 
of  energy  leaving  a wall  segment  j intercepted  by  wall  segment  k. 

number  of  wall  segments 

number  of  fires 


energy  per  unit  time  [W] 

net  heat  flux  leaving  a wall  segment  [W /m^] 

temperature  [K] 


absorbance 


dirac  delta  function;  Sij  = < 


0 

1 


if  i j 
if  i = j 


emittance,  fraction  of  the  black  body  radiation  emitted  by  a gray  sur- 
face or  gray  gas. 


fraction  of  a fire’s  energy  release  rate  that  contributes  to  radiative  heat 
transfer. 


<f>  configuration  factor 

w universal  constant 
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r transmittance 

r average  transmission  factor  where  the  average  is  computed  over  all 

possible  paths  between  two  wall  segments 

cr  Stefan-Boltzmann  constant,  a =■  5.67  x 10~®— 

p reflectance,  fraction  of  the  energy  reflected  by  a surface.  This  fraction 

can  be  related  to  the  emittance,  c , for  a gray  surface  via  e = \ — p. 

T transmittance,  fraction  of  the  energy  passing  through  a gas  unimpeded. 

If  the  gas  has  uniform  absorbency  properties  then  r can  be  computed 
using  the  Beer-Lambert  law  via  where  A is  the  absorbency  of  the 

gas  per  unit  length  and  L is  the  length  of  the  path  through  the  gas. 

LJf^k  a solid  angle;  The  energy  fraction  of  the  f’th  fire  striking  the  k’th  wall 

segment. 


Subscripts 

j,k 

j - k 
f 

f-k 

par 

perp 

total 


wall  segments  j or  k 

from  wall  segment  j to  wall  segment  k 

fire  f 

from  fire  / to  wall  segment  k 

parallel  rectangles  which  are  identical  and  opposite 
perpindicular  rectangles  which  share  a common  edge 
total  energy  release  rate  of  a fire 


Superscripts 

" flux,  a quantity  per  unit  area 

fire  a quantity  due  to  a fire 

in  incoming 

out  outgoing 

L lower  gas  layer 

U upper  gas  layer 

L,gas  a quantity  due  to  the  lower  gas  layer 
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U,  gas 


a quantity  due  to  the  upper  gas  layer 


i,  gas 


a quantity  due  to  the  i’th  gas  layer  where  i can  be  L for  lower  or  U for 
upper  layer 
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A Software  Documentation 


This  section  shows  how  to  use  the  radiation  subroutines:  RAD2,  RAD4  and  RADIO  in  a 
zone  fire  model.  They  model  the  radiation  exchange  between  two,  four  and  ten  wall  segments 
in  an  enclosure  respectively.  These  routines  were  designed  to  be  compatible  with  the  zone 
fire  model  CFAST.  However,  since  these  subroutines  communicate  with  the  calling  routine 
via  the  argument  list  and  not  through  common  blocks  it  should  be  easy  to  incorporate  them 
into  other  zone  fire  models.  They  calculate  the  energy  absorbed  by  the  gas  layers  and  the 
heat  flux  striking  the  wall  segments. 

Linkage  between  the  subroutines  used  by  RAD2,  RAD4  is  represented  by  a computer 
generated  text  file  given  in  Appendix  B.  RAD2,  RAD4  and  RADIO  were  designed  to  be 
modular  so  that  they  can  use  many  of  the  same  routines.  This  simplifies  debugging  and 
shortens  development  time. 

RAD2,  RAD4  and  RADIO  require  wall  properties  such  as  geometric  specifications,  emit- 
tances  and  temperatures;  gas  properties  such  as  absorptivities,  layer  height  and  temper- 
atures. They  calculate  the  rate  of  energy  absorbed  by  the  gas  and  the  radiant  heat  flux 
striking  the  wall  segment. 

DEFFIGS,  LAYFIG,  RMFIG,  PRPFIG2  and  PARFIG2  calculate  the  configuration  fac- 
tors for  RADIO.  RAD2  and  RAD4  use  only  RDPARFIG.  The  subroutine  LAYFIG  and 
RMFIG  implement  the  configuration  factor  formulas  found  in  section  3.1.  LAYFIG  calcu- 
lates the  configuration  factors  between  the  following  10  wall  segments:  ceiling,  4 upper  wall 
segments,  4 lower  wall  segments  and  the  floor.  LAYFIG  is  called  each  time  RADIO  is  called. 
RMFIG  on  the  other  hand  is  called  only  once  per  run.  It  calculates  the  configuration  factors 
between  the  ceiling,  four  walls  and  floor. 

Transmission  factors  are  computed  by  the  subroutine  DEFTAU  and  DEFTAUF.  DEF- 
TAU  computes  transmission  factors  between  wall  surfaces  and  DEFTAUF  computes  trans- 
mission factors  between  wall  surfaces  and  fires.  The  beam  length  chosen  is  based  upon  a 
path  between  the  center  of  two  rectangle. 

The  linear  system  in  equation  (19)  is  solved  directly  or  iteratively  depending  on  how 
close  the  emittances  are  to  1 or  equivalently  how  close  the  wall  surfaces  are  to  a black 
body.  Routines  SGEFA  and  SGESL  are  used  to  solve  the  linear  system  directly.  They  are 
a part  of  Linpack  [17].  For  the  10  wall  case,  a simple  iterative  scheme  based  on  the  Gauss- 
Seidel  algorithm  is  used  to  solve  the  linear  system  when  the  emittances  values  indicate  that 
convergence  can  be  achieved  with  few  iterations.  This  is  not  necessary  for  the  two  and 
four-wall  cases  since  the  speed  up  would  be  inconsequential. 

A.l  Using  Subroutines  RAD4  and  RAD2 

Subroutines  RAD2  and  RAD4  calculate  the  radiation  exchange  between  two  wall  segments 
and  four  wall  segments  respectively.  RAD2  and  RAD4  have  identical  interfaces.  Either 
routine  may  be  used  in  place  of  the  other  without  changing  code  in  the  calling  routine. 
RAD2  returns  the  same  flux  for  the  lower  wall  and  the  floor.  The  following  FORTRAN 
routine  is  an  example  of  how  to  use  RAD2  and  RAD4.  This  routine  was  used  to  generate 
the  data  plotted  in  figure  13. 

PROGRAM  MAIN 
INCLUDE  ’PRECIS. INC’ 

DIMENSION  TWALL(4),  TLAY(2),  EMIS(4),  ABS0RB(2) 

DIMENSION  qFLUX(4),  qLAY(2) 
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c 

c 

c*** 

c 

c 

c*** 

c 

c 

c*** 

c 

c 

c*** 

c 

c 

c 

c 

c*** 

c 


c 

c*** 

c 

11 

c 


c 

c*** 

c 


DIMENSION  QFIRE(IO),  XFIRE(IO),  YFIRE(IO),  ZFIRE(IO) 

OPEN (UNIT=7 , FILE= ’TEST .IN’) 

READ  CEILING,  UPPER  WALL,  LOWER  WALL  AND  FLOOR  TEMPERATURES 
READ(7,*)(TWALL(I),I=1,4) 

READ  UPPER  AND  LOWER  LAYER  TEMPERATURES 
READ(7,*)TLAY(1) ,TLAY(2) 

READ  CEILING,  UPPER  WALL,  LOWER  WALL  AND  FLOOR  EMISIVITIES 

READ(7,*)(EMIS(I) ,1=1,4) 

READ  UPPER  AND  LOWER  GAS  LAYER  ABSORPTIVITIES 

READ(7,*)ABS0RB(1) ,ABS0RB(2) 

READ  ROOM  DIMENSIONS 

READ ( 7 , * ) XROOM , YROOM , ZROOM 

READ  NUMBER  OF  FIRES 

READ(7,*)NFIRE 

DO  11  IFIRE  = 1,  NFIRE 

FOR  EACH  FIRE  READ  FIRE  LOCATION,  AND  FIRE  SIZE 

READ(7,*)XFIRE(IFIRE) ,YFIRE(IFIRE) ,ZFIRE(IFIRE) ,QFIRE(IFIRE) 
CONTINUE 

HLAY  = 3. 

DO  40  I = 1,  101 

COMPARE  RAD2  AND  RAD4  BY  VARYING  UPPER  WALL  TEMPERATURE 

TWALL(2)  = TWALL(3)f(I-l)/100.  + TWALL(1)*(101-I)/100. 

CALL  RAD2( 

I TWALL,TLAY,EMIS, ABSORB, 

I 1, XROOM, YROOM, ZROOM, HLAY, 

I qFIRE,XFIRE,YFIRE,ZFIRE, NFIRE, 

0 QFLUX,QLAY 

) 

WRITE(6,20)TWALL(2) ,(QFLUX(K) ,K=1,4) ,qLAY(l) ,qLAY(2) 

CALL  RAD4( 


I TWALL,TLAY,EMIS, ABSORB, 

I 1, XROOM, YROOM, ZROOM, HLAY, 

I qFIRE,XFIRE,YFIRE,ZFIRE, NFIRE, 
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0 


0 

qFLUX.qiAY 

) 

WRITE(6,20)TWALL(2) ,(qFLUX(K) ,K=1,4) .QLAYd) .qLAY(2) 
WRITE(6,*)’  ’ 

20  F0RMAT(E10.3,1X.9(E10.3,1X)) 

40  CONTINUE 
STOP 
END 

A description  of  the  inputs  follow. 

RAD2,  RAD4  Interface  INPUTS 


TWALL 

TWALL(I)  is  the  temperature  of  the  Fth  surface  [K]  where  1=1, 2, 3, 4 
denotes  the  ceiling,  the  upper  wall,  the  lower  wall  and  the  floor  respec- 
tively 

TLAY 

TLAY(I)  is  the  temperature  of  the  Fth  layer  [K]  where  1=1,2  denotes 
the  upper  and  lower  layers  respectively 

EMIS 

EMlS(l)  is  the  emittance  of  the  ceiling  (1=1),  upper  and  lower  walls 
(1=2,3)  and  floor  (1=4) 

ABSORB 

ABSORB(l)  is  the  absorptance  [m~^]  of  the  upper  (1=1)  and  lower 
layer  (1=2) 

XROOM 

XROOM  is  the  size  of  the  room  [m]  in  the  x’th  coordinate  direction. 

YROOM 

YROOM  is  the  size  of  the  room  [m]  in  the  y’th  coordinate  direction. 

ZROOM 

ZROOM  is  the  size  of  the  room  [m]  in  the  z’th  (vertical)  coordinate 
direction. 

HLAY 

HLAY  is  the  height  of  the  smoke  layer  interface  above  the  floor  [m] 

QFIRE 

QFIRE  is  an  array  of  length  NFIRE;  QFIRE(IFIRE)  is  the  energy 
release  rate  due  to  radiation  of  the  IFIRE ’th  Are  [W] 

XFIRE 

XFIRE  is  the  x coordinate  of  the  fire  location  [m] 

YFIRE 

YFIRE  is  the  y coordinate  of  the  fire  location  [m] 

ZFIRE 

XFIRE  is  the  z (vertical)  coordinate  of  the  fire  location  [m] 

OUTPUTS 


QFLUX 

QFLUX(l)  is  the  radiant  heat  flux  [W/m^]  to  the  I’th  surfaces  where 
i=l,2,3,4  denotes  the  ceiling,  the  upper  wall,  the  lower  wall  and  the 
floor  respectively 
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QLAY 


QLAY(I)  is  the  heat  absorbed  by  the  I’th  layer  where  1=1,2  denotes 
the  upper,  lower  layers  respectively 


A. 2 Using  the  Subroutine  RADIO 

The  interface  for  RADIO  is  essentially  the  same  as  the  interface  for  RAD2  and  RAD4. 
The  variables  TLAY,  EMIS  and  QFLUX  are  arrays  of  size  4 for  RAD2  and  RAD4.  In 
RADIO,  the  arrays  TLAY,  EMIS  and  QFLUX  are  dimensioned  as  size  10  rather  than  size 
4,  corresponding  in  order,  to  the  10  wall  segments:  ceiling,  upper  front,  upper  right,  upper 
back,  upper  left,  lower  front,  lower  right,  lower  back,  lower  left  walls  and  floor.  RADIO 
treats  these  wall  segments  independently. 

B Radiation  Package  Subroutine  Structure  - Who  Calls 
Whom 

This  appendix  contains  information  on  the  subroutines  used  by  RAD2,  RAD4  and  RADIO 
are  related  to  each  other.  This  information  is  provided  to  make  it  easier  to  understand  its 
structure.  Each  entry  has  up  to  four  subheadings:  CALLS,  LIB,  COMMONS  and  CALLED 
BY.  The  subheadings  CALLS  and  LIB  are  similar.  They  both  list  external  references 
to  NAME.,  i.e.,  routines  that  NAME  calls.  The  source  code  for  the  routines  that  are 
listed  under  CALLS  appear  in  the  same  computer  file  as  NAME.  On  the  other  hand,  the 
routines  listed  by  the  LIB  sub-heading  do  not.  Some  examples  of  routines  that  would  appear 
under  the  LIB  heading  are  FORTRAN  supplied  functions  such  as  ABS,  SQRT,  MODE,  etc. 
The  radiation  routines  documented  here  do  not  have  common  blocks  so  there  will  be  no 


information  under  this  heading.  The  routines  that  are  listed  next  to  CALLED  BY  are  those 
routines  that  call  NAME.  This  appendix  was  generated  by  a computer  program  named 
ROADMAP  documented  in  [18]. 

NAME 


CALLS: 


SUBl,  SUB2,  ... 


LIB: 


SUBA,  SUBB,  ... 


COMMONS:  COMl,  COM2,  . . . 

CALLED  BY:  SUBa,  SUBb,  . . . 


SOURCE  ROUTINES 


ROUTINE:  ISAHAX 
LIB: 


ABS 


CALLED  BY:  SGEFA 


ROUTINE:  RAD2 
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CALLS: 

RDABS 

RDFLUX  RDPARFIG  RDSANG 

LIB: 

EXP 

CALLED  BY: 

NONE  - 

NO  ROUTINES  CALL  RAD2 

ROUTINE: 

RAD4 

CALLS: 

RDABS 

RDFANG  RDFLUX  RDFTRAN  RDPARFIG 

RDRTRAN 

I SDOT 

SGEFA 

SGESL 

LIB: 

SQRT 

CALLED  BY: 

NONE  - 

NO  ROUTINES  CALL  RAD4 

ROUTINE: 

RDABS 

CALLED  BY: 

RAD2 

RAD4 

ROUTINE: 

RDFANG 

CALLS: 

RDSANG 

CALLED  BY: 

RAD4 

ROUTINE : 

RDFLUX 

CALLED  BY: 

RAD2 

RAD4 

ROUTINE: 

RDFTRAN 

LIB: 

ABS 

EXP 

CALLED  BY: 

RAD4 

ROUTINE: 

RDPARFIG 

LIB: 

ATAN 

LOG  SQRT 

CALLED  BY: 

RAD2 

RAD4 

ROUTINE: 

RDPRPFIG 

LIB: 

ATAN 

LOG  SQRT 

CALLED  BY: 

NONE  - 

NO  ROUTINES  CALL  RDPRPFIG 

ROUTINE: 

RDRTRAN 

LIB: 

EXP 

CALLED  BY: 

RAD4 

ROUTINE: 

RDSANG 

CALLS : 

RDSANG 1 

LIB: 

ABS 

SIGN 

CALLED  BY: 

RAD2 

RDFANG 

ROUTINE: 

RDSANG 1 

LIB: 

ASIN 

SQRT 

CALLED  BY: 

RDSANG 

ROUTINE: 

SAXPY 

LIB: 

MOD 

CALLED  BY: 

SGEFA 

SGESL 

ROUTINE: 

SDOT 

CALLS: 

SDOT 
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LIB: 

CALLED 

BY: 

HOD 

RAD4 

SDOT 

SGESL 

ROUTINE: 

SGEFA 

CALLS: 

CALLED 

BY: 

ISAMAX 

RAD4 

SAXPY 

SSCAL 

ROUTINE: 

SGESL 

CALLS: 

CALLED 

BY: 

SAXPY 

RAD4 

SDOT 

ROUTINE: 

SNRM2 

LIB: 

CALLED 

BY: 

ABS 

NONE  - 

DBLE  SQRT 

NO  ROUTINES  CALL  SNRM2 

ROUTINE: 

SSCAL 

LIB; 

CALLED 

BY: 

MOD 

SGEFA 

C Subroutine  Glossary 

The  following  is  a glossary  of  subroutines  used  by  the  radiation  packages  RAD2,  RAD4  and 
RADIO. 

BLAS 


LINPACK 


ITER 


LAYFIG 


A collection  basic  linear  algebra  subroutines  used  by  the  Unpack  sub- 
routines SGEFA  and  SGESL  to  solve  the  net  radiation  linear  system 
of  equations.  These  routines  are  documented  in  [19]. 

The  routines  SGEFA  and  SGESL  from  LINPACK  are  used  to  factor 
and  solve  the  net  radiation  linear  system  of  equations.  These  routines 
are  documented  in  [17]. 

This  routine  solves  the  linear  system, 

Ax  = b (28) 

iteratively  by  noting  the  fact  that  A = I — G where 

IIGII  < 1 

= I + G + G^  + G^ + ■■■ 

so  that 

x = iI  + G + G^ + ...)b  (29) 

This  routine  is  only  designed  for  use  with  RADIO.  Note  that  G^  is 
not  explicitly  formed.  Only  a few  iterations  are  required  when  the 
emittances  are  close  to  1. 

This  routine  calculates  the  configuration  factors  for  surfaces  in  a rect- 
angular box  with  a layer  interface.  This  routine  needs  to  be  called  each 
time  the  layer  height  changes. 
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RAD2 

RAD4 

RADIO 

RDABS 

RDFTRAN 

RDRTRAN 

RDFLUX 

RDPARFIG 

RDPRPFIG 

RDSANG 

RDSANGl 

RMFIG 


This  routine  computes  the  radiative  heat  flux  to  the  extended  ceiling 
(ceiling  + upper  wall)  and  the  extended  floor  (floor  + lower  wall)  due 
to  a point  source  Are,  emitting  absorbing  gas  layers  (upper  and  lower) 
and  heat  emitting  wall  segments.  It  also  computes  the  heat  absorbed 
by  the  lower  and  upper  layers. 

This  routine  computes  the  radiative  heat  flux  to  the  ceiling,  upper  wall, 
lower  wall  and  floor  due  to  a point  source  Are,  emitting  absorbing  gas 
layers  (upper  and  lower)  and  heat  emitting  wall  segments.  This  routine 
also  computes  the  heat  absorbed  by  the  lower  and  upper  layers. 

This  routine  computes  the  radiative  heat  flux  to  the  ceiling,  four  upper 
walls  (upper  front,  upper  right,  upper  back  and  upper  left),  four  lower 
walls  (lower  front,  lower  right,  lower  back  and  lower  left),  due  to  a point 
source  fire(s),  emitting/absorbing  gas  layers  (upper  and  lower)  and  heat 
emitting  wall  segments.  This  routine  also  computes  the  heat  absorbed 
by  the  lower  and  upper  layers.  It  performs  the  same  calculations  as 
RAD4  and  RAD2  except  for  the  number  of  wall  segments  considered. 

This  routine  computes  the  energy  absorbed  by  the  upper  and  lower 
layer  due  to  radiation  given  off  by  heat  emitting  rectangles  forming  the 
enclosure.  This  energy  absorption  is  added  to  the  variables,  QLLAY 
and  QULAY  which  were  previously  used  to  contain  the  heat  absorbed 
by  the  lower  and  upper  layers  due  to  gas  emissions  and  fires. 

This  routine  calculates  the  transmissivities,  and  between 

each  fire  f and  each  wall  segment  k. 

This  routine  calculates  the  transmissivities,  and  between  each 
pair  of  wall  segments  j and  k. 

This  routine  calculates  the  ’c’  vector  in  the  net  radiation  equations  of 
Siegel  and  Howell  and  the  enthalpy  absorbed  by  the  lower  and  upper 
layer  fires  due  to  gas  layer  emission  and  fires. 

This  routine  calculates  the  configuration  factor  between  two  identical 
parallel  rectangles. 

This  routine  calculates  the  configuration  factor  between  two  perpen- 
dicular rectangles  that  share  a common  edge. 

This  routine  computes  the  solid  angle  of  an  arbitrary  rectangle. 

This  routine  computes  the  solid  angle  of  a rectangle  whose  corner  is 
tangent  to  the  surface  of  a sphere.  The  center  of  this  sphere  is  the 
location  of  the  point  source  fire. 

This  routine  calculates  the  configuration  factors  for  surfaces  in  a rect- 
angular box.  This  calculation  needs  to  be  done  only  once  per  room. 
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assuming  that  the  room  does  not  change  size  during  the  simulation. 
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